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Abstract 

A search for solar axions has been performed using an axion helioscope which is 
equipped with a 2.3-m long 4T superconducting magnet, a gas container to hold 
dispersion-matching gas, PIN-photodiode X-ray detectors, and a telescope mount 
mechanism to track the sun. A mass region around m a = 1 eV was newly explored. 
From the absence of any evidence, analysis sets a limit on axion-photon coupling 
constant to be <? a77 < 5.6-13.4 x 10 -10 GeV -1 for the axion mass of 0.84 < m a < 
1.00 eV at 95% confidence level. It is the first result to search for the axion in the 
ga-y-y-ma parameter region of the preferred axion models with a magnetic helioscope. 
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1 Introduction 



Quantum chromodynamics (QCD) is the theory of the strong interactions. 
Although QCD has proven remarkably successful, there is a blemish called the 
strong CP problem. The strong CP problem is that the effective Lagrangian of 
QCD has CP violating term but it is not observed, i.e., the experimental value 
of neutron electric dipole moment is smaller than expected by many orders 
of magnitude. Peccei and Quinn proposed an attractive solution to solve this 
problem [T f 2 |5]Hf5] . They introduced a new global U(l) symmetry, Peccei- 
Quinn (PQ) symmetry. When PQ symmetry is spontaneously broken, a new 
effective term arises in QCD Lagrangian which cancels the CP violation term. 
The solution also predicts a new pseudo Nambu-Goldstone boson, axion. The 
expected behavior of an axion is characterized mostly by the scaling factor of 
the PQ symmetry breaking, f a , and so its mass, m a , which is directly related 
to f a by m a = 6 x 10 15 [eV 2 ]// a . 

Axions are expected to be produced in stellar core through their coupling to 
photons with energies of order keV. Especially, the sun can be a powerful 
source of axions and the so-called 'axion helioscope' technique may enable us 
to detect such axions directly [BIT?] . 

The principle of the axion helioscope is illustrated in Fig. [TJ Axions would be 
produced through the Primakoff process in the solar core. The differential flux 
of solar axions at the Earth is approximated by [8 , 9J 



d$ a /d£ = 6.020 x 10 10 [cm- 2 s _1 keV 
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where g a77 is the axion-photon coupling const antPl. Their average energy is 
4.2 keV reflecting actually the core temperature of the sun(~ 3kT), since low 
energy axion production is suppressed due to screening effects[10j. Then, they 
would be coherently converted into X-rays through the inverse process in a 
strong magnetic field at a laboratory. The conversion rate is given by 
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J dzT/2 x J dzB L exp i J dz' (q - y 



(2) 



where z and z' are the coordinate along the incident solar axion, B± is the 
strength of the transverse magnetic field, L is the length of the field along 

1 The formula is different from the one in our previous paper |ll|12j . but both 
coincide numerically with good approximation. 
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z-axis, T is the X-ray absorption coefficient of helium, q = (m 7 — m^)/2E is 
the momentum transfer by the virtual photon, and m 7 is the effective mass of 
the photon which equals zero in vacuum. Eq. d2J) is reduced to 



in case q and B± is constant along z- axis and T = 0. 

In 1997, the first phase measurement [TT] was performed using an axion he- 
lioscope with a dedicated superconducting magnet which is identical to the 
one used in the present experiment, except that the gas container was absent 
and the conversion region was vacuum. Its sensitivity was limited to the ax- 
ion mass region of m a < 0.03 eV due to a loss of coherence by non-zero q in 



If one can adjust m 7 to m a , coherence will be restored for non-zero mass 
axions. This is achieved by filling the conversion region with gas. A photon in 
the X-ray region acquires a positive effective mass in a medium. In light gas, 
such as hydrogen or helium, it is well approximated by 



where a is the fine structure constant, m e is the electron mass, and N e is 
the number density of electrons. We adopted cold helium g clS clS db dispersion- 
matching medium. Here, light gas was preferred since it minimizes self absorp- 
tion by gas. It is worth noting that helium remains at gas state even at 5-6 K, 
the operating temperature of our magnet. Since the bore of the magnet is 
limited in space, the easiest way is to keep the gas at the same temperature as 
the magnet. Moreover, axions as heavy as a few electronvolts can be reached 
with helium gas of only about one atmosphere at this temperature. 

In this way, in 2000, the second phase measurement [12] was performed to 
search for sub-electronvolt axions. This experiment, together with the first 
phase measurement of 1997 [11] with vacuum conversion region, yielded an 
upper limit of g ail < 6.0-10.5 x lO-^GeV- 1 (95% CL) for m a < 0.27 eV. 

In this Letter, we will present the result of the third phase measurement in 
which we scanned the mass region between 0.84 < m a < 1.00 eV using the 
upgraded apparatus to withstand higher pressure gas. 




(3) 



Eq. ©. 




(4) 
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2 Experimental apparatus 



The schematic figure of the axion helioscope is shown in Fig. [21 It is designed 
to track the sun in order to achieve long exposure time. It consists of a su- 
perconducting magnet, X-ray detectors, a gas container, and an altazimuth 
mounting. In the following paragraphs, we will describe each part in due order. 

The superconducting magnet [13] consists of two 2.3-m long race-track shaped 
coils running parallel with a 20-mm wide gap between them. The magnetic 
field in the gap is 4T perpendicular to the helioscope axis. The coils are 
kept at 5-6 K during operation. In order to make it easy to swing this large 
cryogenic apparatus, two devices are engaged. First, the magnet was made 
cryogen-free by making two Gifford-McMahon refrigerators to cool it directly 
by conduction. Second, a persistent current switch was equipped. Thanks to 
this, the magnet can be freed from thick current leads after excitation, and 
the magnetic field is very stable for a long period of time without supplying 
current. 

The container to hold dispersion-matching gas is inserted in the 20 x 92 mm 2 
aperture of the magnet. Its body is made of four 2.3-m long 0.8-mm thick 
stainless- steel square pipes welded side by side to each other. The entire body 
is wrapped with 5N high purity aluminium sheet to achieve high uniformity of 
temperature. The measured thermal conductance between the both ends was 
1 x 10~ 2 W/K at 6 K. One end at the forward side of the container is sealed with 
welded plugs and is suspended firmly by three Kevlar cords, so that thermal 
flow through this end is highly suppressed. Combined thermal conductance of 
three Kevlar cords is estimated to be smaller than that of the container itself 
by four orders of magnitude. Since place-to-place temperature difference of the 
magnet is likely to be much less than IK, the temperature of the container is 
calculated to be uniform along the length within the order of O.lmK, which 
is fairly small and has only negligible effect on the coherence of the axion to 
photon conversion. Radiative thermal flow between the gas container and the 
magnet is the same order of magnitude as the conductive flow through the 
Kevlar cords, and hence negligible, too. The opposite side nearer to the X-ray 
detectors is flanged and fixed to the magnet. At this end of the container, gas is 
separated from vacuum with an X-ray window manufactured by METOREX 
which is transparent to X-ray above 2 keV and can hold gas up to 0.3 MPa at 
liquid helium temperature. 

To have automatic sequential pressure settings of the dispersion-matching gas 
for the scan of the axion mass region around leV, a gas handling system is 
built with 3 Piezo valves (two HORIBASTEC PV1101 and a PV1302) and a 
precision pressure gauge (YOKOGAWA MU101-AH1N). The temperature of 
the gas container was measured by a Lakeshore CGR thermistor. 
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Required pressure for a given mass was determined by the corresponding gas 
density and the temperature based on interpolation of the tables from NIST 
|14j . The error of the target mass is estimated to be less than 5meV by the 
errors of the pressure and the temperature. Since we are scanning a range of 
the axion mass, the error is only crucial to the lower and upper edges of the 
mass range. 

Helium gas is fed to the container through the Piezo valve to have a specified 
pressure setting. If a lower pressure setting is required, then another Piezo 
valve is opened to suck the gas by a vacuum pump connected to the valve. 
Once a proper pressure setting is settled, all the valves are closed and the 
helium gas is kept confined to have a constant electron number density in the 
container until the measurement for the setting is completed. The uniformity 
of the temperature guarantees the homogeneous density along the length of 
the container. The whole process is done step by step automatically to scan 
the axion mass. 

Absorption of X-ray in the helium gas is not negligible and the effect is properly 
calculated in Eq. ([2]). The electron number density might vary slightly along 
the container from the one end to the other because of the gravity when the 
inclination is high, and the coherence could be partly lost accordingly. The 
effect is also taken into account in Eq. (J2]). The decreases of the conversion 
probability P a ^ 7 due to the absorption and the gravity are less than 23% and 
1% , respectively when m 7 is tuned to 1.0 eV at the center of the gas container. 

For emergency exhaust of the gas in case of rapid temperature increase due 
to a magnet quenching, a rupture disk, which is designed to break at 0.248 
MPa, is introduced into the gas handling system to avoid destruction of the 
X-ray window by the over pressure. 

Sixteen PIN photodiodes, Hamamatsu Photonics S3590-06-SPL, are used as 
the X-ray detectors, whose chip sizes are 11 x 11 x 0.5 mm 3 each. In the present 
measurement, however, twelve of them are used for the analysis because four 
went defective through thermal stresses since the measurement of the previ- 
ous phase. The effective area of a photodiode was measured formerly using 
a pencil-beam X-ray source, and found to be larger than 9 x 9 mm 2 . It has 
an inactive surface layer of 0.35 /im [16j . Each chip is mounted on a Kapton 
film bonded to an Invar plate with cryogenic compatible adhesive. The X- 
ray detectors are mounted in a 10-mm thick radiation shielding box made of 
oxygen-free high conductivity copper (OFHC Cu), which is then surrounded 
by a lead shield of about 150 mm thick. The copper shield is operated at about 
60 K, so that it also functions as a cold finger for the X-ray detectors. Details 
on the X-ray detector are given in Refs. [T5][l"6"] . 




The output from each photodiode is fed to a charge sensitive preamplifier 
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whose first-stage FET is at the cryogenic stage near the photodiode chip and 
the preamplifier outputs are digitized using CAMAC flash analog-to-digital 
convertors (FADCs), REPIC RPC-081's, at a sampling rate of 10 MHz. The 
preamplifier outputs are also fed to shaping amplifiers, Clear Pulse CP4026, 
whose outputs are then discriminated to generate triggers. Thus, waveforms of 
the sixteen preamplifier outputs are recorded simultaneously over 50 /j,s before 
and after each trigger to be committed to later off-line analysis. Each detector 
was calibrated by 5.9-keV Mn X-rays from a 55 Fe source installed in front of 
them. The source is manipulated from the outside and is completely retracted 
behind the shield during the axion observations, i.e. during solar tracking. 

The entire axion detector is constructed in a vacuum vessel and the vessel is 
mounted on an altazimuth mount. Its trackable altitude ranges from —28° to 
+28° and its azimuthal direction is designed to be limited only by a limiter 
which prevents the helioscope from endless rotation. However, in the present 
measurement, the azimuthal range is restricted to about 60° because a cable 
handling system for its unmanned operation is not completed yet. The range 
corresponds to an exposure time of about a quarter of a day in observing the 
sun. This is enough for the time being, since background is measured during 
the other three quarters of a day. When the entire cable handling system is 
complete, running time per pressure setting can be shortened by a factor of 
more than two. 

This helioscope mount is driven by two AC servo motors controlled by a com- 
puter (PC). The PC also monitors the azimuthal and altitudinal directions of 
the helioscope regularly by two precision rotary encoders and forms a feed- 
back controlling loop as a whole. The US Naval Observatory Vector Astronomy 
Subroutines (NOVAS) [TFJ were used to calculate the solar position. The alti- 
tudinal origin was determined from a spirit level. While the sun is not directly 
visible from the laboratory in the basement floor, the azimuthal origin was 
first determined by a gyrocompass, which detects the north direction by the 
rotation of the earth within an error of ±8", and then it was introduced to 
the laboratory with a theodolite. 

Since the effective aperture of the helioscope is narrow, it is crucial to de- 
termine its accurate geometry. The axis of the helioscope is defined by two 
cross hairs at the edge of the vacuum vessel. The position of each part of the 
helioscope was measured relative to these cross hairs from their exterior using 
the theodolite when they were installed. The positions of the PIN photodi- 
odes were determined relative to the copper shielding box from a photo image 
taken prior to the installation. As it is hard to estimate analytically the effect 
of the geometrical errors as well as the effect of the size of the axion source, we 
performed a Monte Carlo simulation and found that the overall effective area 
is larger than 371 mm 2 at 99% confidence level. It is smaller than the nominal 
sum of the area of all the living PIN diodes mainly because the line of sight is 
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partly shaded by the walls of the four square pipes of the gas container from 
the axion source region of the sun. Also included are many other small effects 
like geometrical misalignment of the container and the PIN diodes, distortion 
of the square pipes, absorption by thick supporting grids of the X-ray window, 
and so on. 



3 Measurement and Analysis 



From December 2007 through April 2008, a measurement employing dispersion- 
matching gas was performed for 34 photon mass settings with about three days 
of running time per setting to scan around 1 eV, which is shown in Table [TJ 

Before obtaining energy spectra, each event was categorized into two major 
groups, the solar observation and the background. Events while the measured 
direction agreed with the sun are counted as former. When the sun is com- 
pletely out of the magnet aperture, events are counted as latter. Otherwise 
events are discarded. 

We performed numerical pulse shaping to the raw waveforms using the Wiener 
filter. The energy of an X-ray is given by the peak height of a wave after the 
shaping. The shaped waveform is given by 

S*{u)C{u) 

u{u) = "W (5) 



where U(u), S(u), C(u), and N(u) are Fourier transformations of the shaped 
waveform, the ideal signal waveform, the measured waveform, and the noise, 
respectively. Noises are obtained by gathering waveforms while no trigger ex- 
ists, and the ideal signal waveform is approximated by averaging signals from 
5.9-keV X-rays. The response function of this waveform analysis, i.e., non- 
linearity, gain walk by trigger timing, etc., was investigated thoroughly using 
simulated pulses which were obtained by adding the template waveform to the 
noise waveforms. A correction was made based on this numerical simulation. 
Saturation arose at about 25 keV, therefore, E > 20 keV was not used in the 
later analysis. 

Event reduction process is applied in the same way as the second phase mea- 
surement [12] in order to get rid of bad events like microphonics. Applying 
the same cuts to the 55 Fe source data, we found the loss of axion detection 
efficiency by the reduction to be less than 1.5%. 

The background level was about 1.0 x 10~ 5 keV _1 s _1 /PIN(81 mm 2 ) at E = 
5-10 keV. By analysing the calibration data, we found the energy resolution 
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of each photodiode to be 0.8-1.2 keV (FWHM) for 5.9-keV photons. 



In Fig. [31 the energy spectrum of the solar observation with the gas density for 
m 7 = 1.004 eV is shown together with the background spectrum. We searched 
for expected axion signals which scale with for various m a in these spectra. 
The smooth curve in the figure represents an example for the expected axion 
signal where m a = m 7 = 1.004 eV and g ail = 7.7 x 10~ 10 GeV _1 , which 
corresponds to the upper limit estimated as follows. 

A series of least x 2 fittings was performed assuming various m a values. Data 
from the 34 different gas density settings were combined by using the summed 
X 2 of the 34. The energy region of 4-20 keV was used for fitting where the 
efficiency of the trigger system is almost 100% and the FADCs do not sat- 
urate. As a result, no significant excess was seen for any m a , and thus an 
upper limit on g ail at 95% confidence level was given following the Bayesian 
scheme. Fig. [4] shows the limit plotted as a function of m a . Our previous limits 
from the first [TT] and the second [T2] phase measurements and some other 
bounds are also plotted in the same figure. The shown previous limits have 
been updated using newly measured inactive surface layer thickness of the 
PIN photodiode [16]; the difference is, however, marginal. The SOLAX |19j . 
COSME [20] and DAMA [21] are solar axion experiments which exploit the 
coherent conversion (i.e. axion Bragg scattering [52]) on the crystalline planes 
in a germanium and a Nal detector. The experiment by Lazarus et al. [18] and 
CAST [21] are the same kind of experiments as ours. The latter utilizes a large 
decommissioned magnet of the LHC at CERN. Its limit is better than our pre- 
vious limits by a factor of seven in low m a region due to its larger B and L in 
Eq. (J3J). In the region m a > 0.14 eV, however, our previous and present limits 
surpass the limit of CAST. The limit g ail < 2.3 x 10~ 9 GeV _1 is the solar limit 
inferred from the solar age consideration and the limit g ail < 1 x 10 _9 GeV _1 
is a more stringent limit reported by Schlattl et al. |23j based on compari- 
son between the helioseismological sound-speed profile and the standard solar 
evolution models with energy losses by solar axions. Watanabe and Shiba- 
hashi [22] have argued that the helioseismological bound can be lowered to 
g a77 < 4.0 x 10~ 10 GeV _1 if the 'seismic solar model' and the observed solar 
neutrino flux are combined. 



4 Conclusion 

The axion mass around 1 eV has been scanned with an axion helioscope with 
cold helium gas as the dispersion-matching medium in the 4 T x 2.3 m magnetic 
field, but no evidence for solar axions was seen. A new limit on g ail shown in 
Fig. H] was set for 0.84 < m a < 1.00 eV. It is the first result to search for the 
axion in the g a77 -m a parameter region of the preferred axion models [26J with 
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a magnetic helioscope. 

When the complete unmanned operation is ready with the automatic cable 
handling system, the mass scan will be continued to cover still wider mass 
range around leV. 



Acknowledgements 

The authors thank the former director general of KEK, Professor H. Sug- 
awara, for his support in the beginning of the helioscope experiment. This 
research was partially supported by the Japanese Ministry of Education, Sci- 
ence, Sports and Culture, Grant-in-Aid for COE Research and Grant-in-Aid 
for Scientic Research (B), and also by the Matsuo Foundation. 



References 

[I] R.D. Peccei, H.R. Quinn, Phys. Rev. Lett. 38 (1977) 1440. 
[2] R.D. Peccei, H.R. Quinn, Phys. Rev. D 16 (1977) 1791. 

[3] S. Weinberg, Phys. Rev. Lett. 40 (1978) 223. 

[4] F. Wilczek, Phys. Rev. Lett. 40 (1978) 279. 

[5] J.E. Kim, Phys. Rep. 150 (1987) 1. 

[6] P. Sikivie, Phys. Rev. Lett. 51 (1983) 1415; 52 (1984) 695 (erratum). 

[7] K. van Bibber et al., Phys. Rev. D 39 (1989) 2089. 

[8] J. N. Bahcall and M. H. Pinsonneault, Phys. Rev. Lett. 92 (2004) 121301. 

[9] G. G. Raffelt, Contribution to XI International Workshop on "Neutrino 
Telescope", Venice, Italy, Feb. 22-25. 2005. 

[10] G. G. Raffelt, Phys. Rev. D 33 (1986) 897. 

[II] S. Moriyama et al., Phys. Lett. B 434 (1998) 147. 
[12] Y. Inoue et al., Phys. Lett. B 536 (2002) 18. 

[13] Y. Sato et al., Development of a Cryogen-free Superconducting Dipole Magnet, 
in: Proc. of the 15th International Conference on Magnet Technology (MT-15) 
(Beijing, October 1997), eds.: L. Liangzhen, S. Guoliao, Y. Luguang (Science 
Press, Beijing, 1998) pp. 262-265; KEK-Preprint-97-202 (November, 1997). 



9 



[14] Vincent D. Arp and Robert D. McCarty, Thermophysical Properties of Helium- 
4 from 0.8 to 1500 K with Pressures to 2000 MPa, NIST Technical Note 1334, 
(U.S. Department of Commerce, National Technical Information Service, 1989). 

[15] T. Namba et al., Nucl. Instr. Meth. A 489 (2002) 224. 

[16] Y. Akimoto et al., Nucl. Instr. Meth. A 557 (2006) 684. 

[17] G. H. Kaplan et al., Astronomical Journal 97 (1989) 1197; 
URL: |http : //aa . usno . navy . mil/A A/ software/ novas/. 

[18] D. M. Lazarus et al., Phys. Rev. Lett. 69 (1992) 2333. 

[19] A. O. Gattone et al., Nucl. Phys. B(Proc. Suppl.)70 (1999) 59. 

[20] A. Morales et al., Astropart. Phys. 16 (2002) 325. 

[21] R. Bernabei et al., Phys. Lett. B 515 (2001) 6. 

[22] E. A. Paschos, K. Zioutas, Phys. Lett. B323 (1994) 367. 

[23] H. Schlattl, A. Weiss, G. Raffelt, Astropart. Phys. 10 (1999) 353. 

[24] K. Zioutas et al., Phys. Rev. Lett. 94 (2005) 121301, C. Eleftheriadis, 
arXiv:0706.0637vl, and S. Andriamonje et al. (CAST Collaboration), 
J. Cosmol. Astropart. Phys. 04 (2007) 010. 

[25] S. Watanabe, H. Shibahashi, hep-ph/01 12012 

[26] D. B. Kaplan, Nucl. Phys. B 260 (1985) 215, M. Srednicki, Nucl. Phys. B 260 
(1985) 689, and S. L. Cheng, C. Q. Geng, W. -T. Ni, Phys. Rev. D 52 (1995) 
3132. 



10 



In the solar core 



In the magnet 




Fig. 1. The solar axions produced via the Primakoff process in the solar core are, 
then, converted into X-rays via the reverse process in the magnet. 
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Fig. 3. The left figure shows the energy spectrum of the solar observation (error bars) 
and the background spectrum (solid line) for the effective PIN photodiode area of 
371 mm 2 when the gas density was tuned to m 7 = 1.004 eV. The right figure shows 
the net energy spectrum of the left where the background is subtracted from the 
solar observation. The solid line shows the expected solar axion energy spectrum. 
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Fig. 4. The left figure is the exclusion plot on g aTy to m a . The new limit and the pre- 
vious ones |ll|12] are plotted in solid lines. Dashed lines are the limit by Lazarus et 
al. [18], the limit by CAST experiment [24] . the limit by SOLAX experiment |19j . 
the limit by COSME experiment [20], the limit by DAMA experiment [21], the 
limit inferred from the solar age consideration, and the helioseismological bound. 
The hatched area corresponds to the preferred axion models [26J. The right figure 
shows the magnified view of the new limit. 
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oou.y 


n q^7 
u.yo i 


zUbzb 


booUz 


00 ( .u 


n qr9 
u.yoz 


99/I 78 

ZZ4 / 8 


7QA/J 1 

i oUbl 


^R9 n 

OOZ.U 


u.yoo 


991 Ofi 
ZZIZO 


bboy4 


00 1 .0 


n Q71 
u.y 1 1 


OA 7on 
z4 1 yU 


7/1 c c a 
/ 4000 


^71 1 
Olo.l 


n Q7R 

u.y i u 


9zLzLfizL 


uyuyu 


578 5 


Q80 


23120 


69297 


584.0 


0.985 


23922 


73955 


588.4 


0.988 


29533 


91299 


597.0 


0.996 


35946 


110559 


602.6 


1.000 


36503 


108174 


607.3 


1.004 


32584 


97671 



Table 1 

Table of the gas settings and each live time. 
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